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We present a method to calculate the gravitational energy when asymptotic bound- 
JH , ary conditions for the space-time are not given. This is the situation for most of the 

cosmological models. The expression for the gravitational energy is obtained in the con- 
text of the teleparallel equivalent of general relativity. We apply our method first to 
the Schwarzschild-de Sitter solution of Einstein's equation, and then to the Robertson- 
Walker Universe. We show that in the first case our method leads to an average energy 
density of the vacuum space-time, and in latter case the energy vanishes in the case of 
null curvature. 

Keywords: Teleparallelism; Torsion Tensor; Cosmology. 

1. Introduction 

Recently some investigations have revived the question about the total energy of 
the Universe, including the material energy and the energy of the gravitational field 
LUlZlloJ|4|lo|lDjL£J This question is nontrivial because the usual treatment of the energy 
of a cosmological model with pseudotensors or Komar's integral does not seem 
appropriate, since the model is described by a metric that is not asymptotically flat 
^1 The problem holds even in an empty Universe described by de Sitter's solution, 
since in this case there is a (non-flat) cosmological horizon. 
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2 Ulhoa et. al 

The motivation for finding an expression for the energy of the universe comes 
from the fact that this subject may help to answer some questions about the origin 
and evolution of the universe, with consequent improvements for the cosmological 
models ^1. 

Several investigations and results indicate that the total energy of the universe is 
zero. However these results are dubious since most of them are obtained by means 
of pseudotensors, which are coordinate-dependent expressions whose geometrical 
meaning is not clear. The results obtained by means of Komar's integrals, on the 
other hand, depend on the normalization of a Killing vector and in the case of cos- 
mological systems, which are not asymptotically flat, there is no physically preferred 
choice for such vectors. 

In the framework of the Teleparallel Equivalent of General Relativity (TEGR) 
the gravitational energy is well defined for finite volumes of the three dimen- 
sional space, and in particular for asymptotically flat space-times. When asymptotic 
boundary condition are not available, for instance in the case of the Robertson- 
Walker model and of the de Sitter Universe, a new method has to be used. As 
remarked by Faddeev an expression for the gravitational energy must vanish 
only in the total absence of matter and gravitational field. Therefore in order to 
avoid the use of boundary conditions , it is possible, and perhaps mandatory, to use 
regularized expressions which do not require the establishment of asymptotic 
conditions. The purpose of this paper is to point out that the use of a regularized 
expression for the gravitational energy is quite suitable to cosmological models. 

This paper is organized as follows. In section 2, we introduce the formalism of 
teleparallel gravity. We show how to define the gravitational energy in this context 
and present the definition of the regularized expression for the gravitational energy, 
which can be applied for arbitrary field configurations. In section 3, we propose a 
method to evaluate the energy based on the regularized expression, since in cosmol- 
ogy the lack of spatial asymptotic conditions prevents the application of ordinary 
methods. We evaluate the energy for two configurations which have relevance in 
cosmology, the Schwarzschild-de Sitter's solution and the Robertson- Walker model 
for the universe. Finally we present some concluding remarks. 

Notation: space-time indices /x, z/, ... and SO(3,l) indices a, 6, ... run from to 
3. Time and space indices are indicated according to fi = 0,i, a = (0), (i). The 
tetrad field is denoted by e a M , and the flat, Minkowski space-time metric tensor 
raises and lowers tetrad indices and is fixed by r\ a b = e^e^g^ = ( — 1, +1, +1, +1). 
The determinant of the tetrad field is represented by e = det(e a ^) and we use the 
constants G = c = 1. 

2. The regularized gravitational energy-momentum expression in 
TEGR 

The TEGR is an alternative formulation of the theory of general relativity, and 
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is constructed out of tetrad fields in the Weitzenbock space-time of distant paral- 
lelism. In the Hamiltonian formulation the constraint equations are interpreted as 
energy, momentum and angular momentum equations for the gravitational field 
Two vectors are said to be parallel if their projections on the tangent space by 
means of the action of the tetrad field are equal. Thus, considering two vectors, 
V a {x) = e a u V»(x) and V a (x + dx) = e a a V^{x) + (e a x d^V x + V x d^e a x)dx^ = 
V a + e a ^(V \V^)dx x , distant from each other by an infinitesimal displacement, the 
teleparallelism or distant parallelism is obtained if 



V U V» = dyV 1 + {e a »d v e aX )V x = , 

where e afJ, d u e a \ — r^' A is called the Weitzenbock connection. It is simple to check 
that V^e" x = 0. 

It is well known that in the Riemannian geometry the Christoffel symbols °r* 
are symmetric in the lower indices and the corresponding torsion tensor vanishes. 
However, the field equations in teleparallel gravity are constructed out of the torsion 
tensor T x M „ of the Weitzenbock connection, where T x ^ = e a x T a and 



T\v(e) = d ti e a v -d u e%. (1) 

The curvature tensor constructed out of the Weitzenbock connection = 
e aX dfj,e aL , vanishes identically. 

Let °u>f la b represent the torsion-free Levi-Civita connection, 

a ^->fj,ab = — 7J eC fJ,(^abc — Qbac — ^cab) , (2) 



^abc — &av{pb ^ &fj,&c C c ^ O^Cft ) . 

3 °r^ J/ and the Levi 
by d M e a „ + V%e 6 „- °T x e a x = 0, or 



The Christoffel symbols °r* and the Levi-Civita connection are identically related 



\ fJiV = e o^eav + e ( Uuabje v . 

In the teleparallel equivalent of GR the following identity is relevant in the con- 
struction of the Lagrangian density, 

°^uab = —K^ab , (3) 

where K^ ab = \e a x e b V (T\ MV + T v \ a + T^v) is the contorsion tensor. 

Let us first write the curvature scalar constructed out of ([3]). It is possible to 
show that such a quantity is given by 

eR(°u) = -e (-T abc T abc + ^T abc T bac - T a T a j + 2d„(eT») . (4) 
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where e is the determinant of the tetrad field e a M and T a — T b b a . Both sides of 
(U) are invariant under Lorentz transformations. Dropping the divergence term we 
construct the Lagrangian density 

C(e ap ) = -ke (^T abc T abc + ^T abc T bac - T a T a - 2A^) - C M 

= -k eZ abc T abc - Cm + 2keA , (5) 

where k = l/(167r), Cm stands for the Lagrangian density for the matter fields, A 
is the cosmological constant and T, abc is defined by 

Y^abc ^_^rjiabc _|_ rpbac rjicab^ _|_ ^_ ^q acr J 1 ^ Tj ab T C ^j (6) 

Performing the variational derivative of the Lagrangian density with respect to 
the tetrad field e a \ we get, after some rearrangements, the field equation 

d v {e£ aXv ) = -^e e a + T Xfl ) + \^ee aX , (7) 

where 

= k(AT, bcX T bc » - g x ^ bcd T bcd ) , (8) 

and e a p T Aai = \^r^ corresponds to the energy-momentum tensor of matter fields. 
It is possible to show, by explicit calculations, that ([7]) is equivalent to the Einstein's 
field equations with the cosmological term. 

Now let us analyze the meaning of t x ^. In view of the ant-symmetry property 

£ap V = j t f U ows that 



dx [ee\(t Xfi +T ,Xfi )] =0, (9) 

which is a local balance equation, where T lX ^ = T x ^+2kkg x ^ . Therefore we identify 
t XfJ " as the gravitational energy-momentum tensor 021. 

The integration of t x » + T' x » over a hyper-surface x° = constant defines the 
energy-momentum vector due to gravitational and matter fields, 

P a = [ d 3 xee a fl {t°^+T' ^), (10) 
Jv 

where V is a volume of the three-dimensional space. Since the integrand is a scalar 
density under coordinate transformations, V may be finite as well as infinite (in 
which case S is finite or S — > oo, respectively). In view of the field equations, eq. 
(fTU| may be written as 
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P a = - I d 3 xdjll aj = - f dSj n aj , (11) 
Jv Js 

where IP- 7 = — 4/ce T, a0: > . It is interesting to note that the above expression is invari- 
ant under coordinate transformations of the three-dimensional space, under time 
reparametrizations, and is projected on the tangent space which means that it is 
coordinate independent. We point out that the definition above was first presented 
inRef. DJ. 

In this formalism gravitation is considered as a manifestation of torsion. There- 
fore in the flat space-time the torsi on t ensor is supposed to vanish. However this con- 
dition in some cases is not satisfied and in this case we have to adopt regularized 
expressions. A regularized expression for the gravitational energy-momentum is 
defined by 



Prea = ~ I d 3 x d k [U ak (e) - IP* (E)} , (12) 
Jv 

where LP- 7 (E) is constructed out of flat tetrads E a M which are obtained from e a M 
by requiring the vanishing of the physical parameters. This definition guarantees 
that the energy-momentum of the flat space-time always vanishes. In the context 
of the energy problem in general relativity, the i dea an d concept of a regularized 
expression was first considered by Brown and York H3H5] m rea lm of the quasilo- 
cal definition of the gravitational energy. Conceptually, our approach regarding the 
regularization does not differ from that considered by these authors. 

In eq. (fT2|) e a M and E a M are separately solutions of Einstein's equations. The 
flat tetrads E a ^ may yield nonvanishing values for the tensor T° and for this 
reason the regularization procedure is necessary. Note that definition (|10j) follows 
from Einstein's equations, and a similar equation holds for E a M . Thus eq. (fT2|) may 
be understood as the subtraction of the two Einstein's equations for e a M and E a ^ . 
We remark, however, that in view of eq. (j4|) a local SO(3,l) transformation of the 
tetrad field yields only a surface term in the Lagrangian density ([5]). The surface 
term does not contribute to the field equations, which are unaffected by a local 
SO(3,l) transformation of the tetrad fields. 

Tetrad fields that in the flat space-time limit reduce to unit vectors in the r, 
9 and <p directions yield nonvanishing torsion tensor components The great 
advantage of using such tetrad fields is that they do not require the prescription of 
boundary conditions. This is precisely the case in cosmology. However, they require 
the use of a regularization procedure. 



3. The Procedure to Evaluate the Gravitational Energy 

For any solution of Einstein's field equations, such as the Schwarzschild or Kerr 
solutions, the choice of the tetrad field fixes the frame adapted to a certain observer 



18, 2010 0:19 WSPC/INSTRUCTION FILE GEPforCM 



6 Ulhoa et. al 



in space-time. Such a choice may be motivated, for instance, by the asymptotic 
conditions of the tetrad at spatial infinity, in which case we take the tetrad field to 
represent the flat space-time in this limit. However, in cosmology this is not the usual 
situation, since one does not require boundary conditions for the metric tensor as one 
does in the analysis of localized material systems in gravitation. Instead of boundary 
conditions, in cosmology one usually requires the concept of homogeneity of the 
three-dimensional space, in agreement with the cosmological principle, which asserts 
that we do not occupy any special position in the universe. The concept of isotropy 
of space is an additional requirement of the cosmological models. Thus in cosmology 
we consider the tetrad field to represent a local basis for an arbitrary observer in 
the chosen coordinate system, and make use of the regularization procedure. Let us 
analyze the procedure in the following cases. 



3.1. The Schwarzschild-de Sitter Space-Time 

The most general spherically symmetric vacuum solution of the field equations 
with a positive cosmological constant A is the Schwarzschild-de Sitter solution, 



ds 2 = - - ^ - ^) dt 2 + (l-^-^j 1 dr 2 +r 2 d6 2 +r 2 sin 2 Odtf , (13) 

where R = \J~^- This metric specifies the de Sitter vacuum solution when we set 
m = 0, and the Schwarzschild solution when A = 0. 
Let us choose a diagonal tetrad field like 



e%(t,r,M) 



( A 
A' 1 
r 

\0 rsin(9/ 



(14) 



where A 



1 



1/2 



The determinant e of ([14]) reads e 



r sin^. In 

=rd6 



terms of one-forms the spatial components are given by e^ 1 -* = A~ x dr 
and e*- 3 -* = r sin 9 dcf>. The spatial components of the tetrad field represent unit 
vectors in the r, 6 and <p directions. Therefore such a tetrad field is interpreted as 
the usual basis in spherical coordinates for local observers. 

In order to calculate the energy of the configuration defined by (fl3|) , we have to 
evaluate the component S 001 , obtained from ([5]), which after algebraic manipula- 
tions is given by 



^001 



1 



9 9 \9 J212+5 1 3i3) 



(15) 



where the components of torsion tensor appearing in the above expression are 
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T212 = r , 

T313 = r sin 2 6 . (16) 
Thus integrating over a surface of constant radius ro in pip it follows that 



P(°) = 4k j> cfSi(e£ 001 ) 



is 

= lim (-rV-5oo) ■ (17) 



However, using the regularized expression given by eq. (|12|) . we find that the energy 
contained within a surface of constant radius ro is 




(18) 

Let us evaluate expression (fl8|) for the range of values of ro such that 



2m r 2 
— « X ' 152 «! 



The above conditions mean that the cosmological constant is very small. Ex- 

^ and ^ 



panding (|T5|) and neglecting all powers of both — and ^j- we arrive at 



P^ g =m+ r fA. (19) 

which yields, for A = 0, the well known value of the energy of the Schwarzschild 
space-time. 

Now we will instead consider the vacuum de Sitter solution only, since in this 
case the main features are not altered by the introduction of a mass m at r = 0. 
Let us assign Eds as the value of the energy in the absence of the mass m, which is 
the background (vacuum) energy. 

The total gravitational energy E,is contained in the physical region of the vac- 
uum de Sitter space is obtained from (fT8]) by making m — and r = R. It follows 
that 



E dS =R=^j. (20) 

If we substitute r = Rsm\ in (fT5|) . then we shall have a space where the coordinate 
X varies from to 7r. The surface of a sphere in these coordinates is S — Ai:R 2 sin 2 x, 
which increase reaching its maximum value 4-kR 2 at \ — 7r/2, after which it decrease 
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to zero at \ = 7rE3. The volume is V = f£ dx^irR 3 sin 2 x therefore the mean energy 
density is given by 

Eds _ E d s A 
V ~ 27T 2 i? 3 " 6tt 2 ■ [ ' 

Thus when we consider an Universe filled with matter, the accelerated rate of ex- 
pansion of the Universe, which is usually explained by the hypothesis of the so 
called dark energy, may be due to the mean energy density of the vacuum space- 
time. Expression ([21]) already appeared in Ref. [I7J however we decide to reconsider 
the same problem since in the present procedure we do not specify any boundary 
conditions for the tetrad field. This method is better suited to describe cosmological 
models. 



3.2. The Homogeneous and Isotropic Universe 

The cosmological principle asserts that the large-scale structure of the Universe 
reveals homogeneity and isotropy^l The most general form of a line element that 
preserves such features may be written as^[ 



ds 2 = -dt 2 + a 2 (t) 



T ^ + ,W + sinW) 



(22) 



where a(t) = S(t)/\K\^ 2 if K ^ and a(t) = S(t) if K = 0. S(t) is the scale factor 
and K is the constant curvature of space. Here K = \K\k' where k' assumes the 
values +1, 0, —1 which correspond to (i) a three-space of constant positive curvature, 
(ii) a flat space or (iii) a space of constant negative curvature, respectively. By 
introducing a new radial parameter f related to r by 



(23) 



l + k'f 2 /4 
we can write the line element above as 

ds " = - dt2 + T, ^K,^ ( df2 + f2 (de 2 + sin 2 0# 2 )) . (24) 
(1 + k'f 2 /4) 2 K v " v ' 

This second form is called the Roberston- Walker (RW) line element. In the following, 
we will consider the first line element. The quantity a(t) is obtained as a solution 
of the Friedman's equations, 

a 2 k! 

3^ + 3^ = 87rp, (25) 
a a 

d 2 „ a kl „ 

— +2— + — = -8irp, 26 
a z a 1 a z 

where p and p denote the density of the cosmological medium and its pressure. 
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In what follows we will show that in the context of the TEGR it is possible to 
arrive at an expression for the gravitational energy of the universe for closed, open 
and flat cases, respectively, using the following diagonal tetrad field, 



<(*,r,M) = 



OA 





ra(t) 
\0 ra(t) sin0 J 



(27) 



°(«) 



The determinant e of (|2~T1) reads e = Aa 2 (t)r 2 sin 9. 



where A =- 

In order to obtain the energy of this system, we have to calculate E 001 , which 
reads 



s 001 = ^00 5 ll (5 22 T2i2 +5 33 T3i3)j (2g) 

where 



T212 = ra 2 (t) , 

T 313 = ra 2 {t) sin 2 9. (29) 

Taking into account eq. (|11[) . we integrate now expression (1281) over a two- 
dimensional spacelike surface S defined by a radius r = vq = constant. We arrive 
at, 

P(°) = Akj dSi(e£ 001 ) 

= -aroVl-fc'M 2 ■ (30) 
We apply now eq. (fT2"j) , and after the regularization procedure we obtain 



= aro (l - y/l ~ fc'M 2 ) . (31) 

The expression (|3ip was obtained in Ref. EElby means of a totally different method. 
This fact supports our own procedure and result. 

In the case of flat space-time k' = 0, the energy is zero. This result is consistent 
with the fact that in this case we have an open (infinit) space with null curvature. 
This means that in this configuration, the gravitational energy is minus the energy 
of the fields of matter for any value of volume considered. In this configuration, 
the gravitational energy and the energy of matter fields exactly cancel out and the 
resultant space is flat. 

For the case of k' — 1 we have a closed space with constant positive curvature, 
and the total energy given by eq. (|31[) for an arbitrary volume of space is always 
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positive and linear in the scale factor. This result reveals the fact that, in this 
case, the gravitational energy and the energy of the matter fields do not exactly 
cancel out, and the difference is likely to be responsible for the nonvanishing positive 
curvature. 

Finally in the case of space with constant negative curvature k' = — 1 (open 
space) we can see from eq. (|3Tj) that the resultant energy (gravitational energy plus 
the energy of matter fields) is always negative for any value of volume, in particular 
it is infinite negative when ro — ¥ oo. Here it is important to note that in this case 
we have an infinite space with constant negative curvature. It is likely that in this 
case, the total negative energy, is responsible for the constant negative curvature of 
the space. 

Recently ^3 negative energy was obtained for cosmological models describing 
open universes. In [20] all 9 types of Bianchi homogeneous cosmological models has 
been examined by using an energy expression obtained out of a Hamiltonian bound- 
ary term in the context of tetrad-teleparallel gravity. In contrast to our results, the 
authors have found zero energy for some non-flat cosmological models. However for 
space with negative curvature the results in [20] have the same sign of our result for 
space with negative curvature. As stated in [20] for the case k' = — 1, each region of 
constant negative curvature acts like a concave lens exactly as if it had a negative 
matter density repelling light rays from these. 

4. Conclusion 

It is well known that solutions of Einstein's equations that describe cosmological 
models do not have well defined boundary conditions in the same way as the config- 
urations of isolated sources. The absence of such conditions in these configurations 
is a difficulty for the determination of a definite expression of gravitational energy 
for cosmological models. In this work we have proposed a coordinate independent 
method to investigate the energy of gravitational field configurations when asymp- 
totically flat conditions are not available. The definition given by eq. (fP2|) in the 
absence of matter fields gives the gravitational energy, and in the presence of matter 
fields gives the total energy. 

We applied our method to the homogeneous and isotropic Universe and to the 
Schwarzschild-de Sitter solution. When the cosmological constant is set to zero in 
the latter case, it is well known that the space-time is asymptotically flat. When 
the cosmological constant is nonvanishing, the cosmological horizon prevents the 
existence of such boundary (asymptotic) conditions. When one tries to use the 
tetrad field as a local frame such that the spatial components are unit vectors along 
the r, 9 and <f> directions, which is realized by means of a diagonal tetrad of the 
type given by eq. (|14j) . some drawbacks may arise. For example, the emergence of 
non-vanishing torsion components in Minkowski space-time. These drawbacks are 
eliminated by means of the r egul arization procedure, which is similar to that given in 
the work by Brown and YorkL^lfor quasi- local energy expression. The regularization 
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procedure is powerful because it dispenses the use of boundary conditions for the 
gravitational field. 
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